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Abstract 
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A brief derivation of a specific regularization for the magnetic gas dynamic system of 
| 7 | i equations is given in the case of general equations of gas state (in presence of a body force 

iy-j " and a heat source). The entropy balance equation in two forms is also derived for the 

system. For a constant regularization parameter and under a standard condition on the 
heat source, we show that the entropy production rate is nonnegative. 
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1 Introduction 



Regularized (or quasi-) gas dynamic systems of equations are exploited for 
various purposes including the construction a class of so-called kinetically con- 
sistent finite-difference methods for gas dynamics simulations. The correspond- 
ing background, designing of the finite-difference methods and various applica- 
tions are presented in monographs |CE]-[3]- 

Several new issues in the mathematical treatment of this approach have been 
recently developed in In particular, the case of a general equations of 

gas state (in presence of a body force and a heat source) have been covered 
in [7]- [10] where the law of non-decreasing entropy, the Petrovskii parabolicity 
and the linearized stability of equilibrium solutions have been established. 

In Magneto Gas Dynamics (MGD) [T2i 113], an application of this approach 
has been recently given. In [TdT [TB] the corresponding regularized system has 
been formally written in the case of a perfect polytropic gas (in absence of 
body forces and heat sources) and some successful numerical results has been 
presented in ID and 2D cases aimed at astrophysical applications. 
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In this paper, following a recent formalism from [TT], we give a brief complete 
derivation of the regularized MGD system of equations in the standard form of 
mass, momentum and total energy balance equations together with the Faraday 
equation covering the case of general gas state equations (in presence of body 
force and heat source). Our formulation of the equations is more standard and 
seems more suitable for further discretization. In addition, we present a useful 
regularized internal energy balance equation. 

It is well known that a crucial point in the physical and mathematical (see 
14] for a complete account of this last point) correctness of a gas dynamics 



system is an adequate entropy balance equation, and our main result is the 
derivation of such an entropy balance equation for the regularized MGD system 
of equations. We write it down in two forms enlarging the recent corresponding 
results from j7], P, [TO], moreover, for a constant regularization parameter and 
under a standard condition on the heat source, we prove that the corresponding 
entropy production rate is nonnegative. 

2 A regularization of the magnetic gas dynamics system 
of equations 

We begin with the classical Navier-Stokes system of equations for a viscous 
compressible gas flow taking into account the magnetic field, a body force and a 
heat source. The system consists of the mass, the impulse and the total energy 
balance equation together with the Faraday equation and the equation of the 
absence of magnetic charge 

dtp + div(/m) = 0, (1) 
dt(pu) + div(pu (g)u-B(g)B) + V(p+ ||B| 2 ) = dWU NS + pF, (2) 
d t (E + ±|B| 2 ) + div [{E + p+ |B| 2 )u - (u ■ B)B] 

= div(xW + U NS u) + pu ■ F + Q, (3) 
d t B + div(u ® B - B <8> u) = 0, (4) 
divB = 0. (5) 

We consider the gas density p > 0, the velocity u = (ui, . . . ,u n ), the ab- 
solute temperature 6 > and the magnetic field strength B as the basic un- 
known functions (the multiplier -|= is included into the definition of B). In 

1 |2 

addition, the equations include the total non-magnetic energy E = - p\u\ + pe, 
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the pressure p and the specific internal energy e. The system is considered for 
(x, t) G Q x (0, T), where Q is a domain in R n , n ^ 1. 

Concerning the notation, hereafter the operators div and V = (di, . . . , <9 n ) 
are taken with respect to the spatial variables x = (xi, . . . ,x n ). Also di and 
dt are the partial derivatives in Xi and t. The divergence of a tensor is taken 
with respect to its first index. The signs (g) and • denote the tensor and inner 
products of vectors, and in the inner products such as u • the sign • is 
omitted for brevity. Also the sign : means the inner product of tensors. 

We take general state equations 

P = p(p,9), £ = e(p,0) (6) 
linked by the Maxwell relation 

p = 6p e + p 2 e p (7) 
and satisfying the thermodynamic stability conditions 

p p ^ 0, e e > 0. 

Hereafter p p , p@, e p and eg are partial derivatives of the state functions (Ej). 

In the above equations, Uns is the classical Navier-Stokes viscous stress 
tensor 

U NS = n NS (u) = /i [2B(u) - § (divu)l] + A(divu)I, D(u) = \ (Vu + Vu T ) 

with the dynamic viscosity coefficient \i = fi(p, 0) ^ and the bulk viscosity 
coefficient A = X(p, 9) ^ 0, also I is the identity tensor (of order n) and Vu = 
{d{Uj}fj =1 . Moreover, x = x(p,6) ^ is the heat conductivity coefficient, the 
given functions F = F(x, t) and Q = Q(x, t) ^ are the density of body forces 
and the power of heat sources. The magnetic viscosity is neglected. 

To regularize this system of equations, in general we follow [21 E] (see also 
| I17| ) but apply the very recent simpler formalism from [TT] and replace terms 

pa, pu (g) u - B (2) B, V (p+ ||B| 2 ) - pF, 

(E + p+ |B| 2 )u- (u- B)B, pu-F, u<2)B-B(8)u 

in the divergent summands in equations (00), (EJ), (EJ) and (Hj) respectively by 

pu + rdt(pu), pu0u — B(g)B + rd t (pu <8) u — B B), 

V[p+l\B\ 2 + Td t (p+±\B\ 2 )} - {p + rd t p)F : 

(E + p+ |B| 2 )u - (u • B)B + rd t [{E + p+ |B| 2 )u - (u • B)B], 

[pu + rdt(pu)} -F, u®B-B(g)u + rd t (u ® B - B <g) u) 
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with a relaxation parameter r = r(p, e, u, B) > 0. Here the hat : over the 
derivative d t means that it is calculated by virtue of the equations neglecting 
viscosity and heat conductivity (i.e., for zero Uns and x) as for the Euler MGD 
system of equations. 

Proposition 1. The regularized MGD system of equations has the form 

d t p + div[p(u - w)] = 0, (8) 

d t (pu) + div[p(u -w)®u-B(g)B] + V(p+ ||B| 2 ) 

= div n + [p - t div(pu)]F, (9) 

d t (E + i|B| 2 ) + div [(E + p)(u - w) + |B| 2 (u - w) - ((u - w) • B)B] 

= div [-q + t((3 ■ B)u + nu] + p(u - w) • F + Q, (10) 

d t B + div[(u - w) <g> B - B <8> (u - w)] = div[r(u <8> (3 - (3 <8> u)] , (11) 

divB = (12) 

provided that divB| i= o = 0. 

Here the auxiliary velocity vector- functions w and w are given by formulas 

T 



w = - [div(/m <g) u - B ® B) + V (p + ±\B\ 2 ) - pF], 
p 

w = - [p(uV)u - div(B <g> B) + V (p + ||B| 2 ) - pF}. 
P 

The non-symmetric regularized viscous stress tensor has the form 

n = U NS + pu <8> w - t(/3 &> B + B <g> /3) 



(13) 
(14) 



+r ^uVp + /?C 2 div u + /3 • B 



pee 



Q I, 



(15) 



where C s ^ is the speed of sound in the gas defined by a known formula 

free 



(16) 



(for example see [T8^). The regularized heat flux q is given by a known formula 

V 



p uVe 5 uVp - Q 



P^ 



- q = kS79 + r 

and the auxiliary vector-function P has the form 

P = div(u(g)B - B <8>u). 



u. 



(17) 



(18) 
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Proof. By virtue of equations (pQ) and (j2j) we have 

d t p = -div(pu), (19) 

r4(pu) = -pw (20) 

with w given by ( 1131) . The latter formula implies the regularized mass balance 
equation ((HI). Also the following formula 

rd t (pu u) = T{d t (pu) Cg u + u Cg [4(pu) - (4p)u]} 
= -pw Cg u — u Cg [pw — r div(pu)u] = — pw (g) u — u Cg pw (21) 

holds taking into account that 

pw = pw + t div(pu)u (22) 

with w given by (TT4l). Notice also the related useful formula 

rdtu = - d t (pu) - - t p)u = -w, (23) 
p p 

see (EUj) and (Hi). 

The well known equation 

d t p = -( uVp + pC 2 s div u - — Q ] (24) 
V P £ e J 

holds, where the speed of sound C s is given by (TTBj). The right-hand side of the 
equation does not depend on B. 
Since clearly 

4B = -(3 (25) 

with (3 given by (EES), we get 

4(B(g)B) = -(^(g)B + B(g/3), 4Q|B| 2 ) = -/3-B, (26) 

and after recalling (E3j) 

r4(u (g) B - B <g) u) = -[w(gB-B(gw + r(u(g/3-/3(gu)]. (27) 

Collecting the above formulas (ED), (Elj) and (ESJ) , we derive the regularized 
impulse balance equation ([9]) with the regularized viscosity tensor in the form 
(H5j). Also formula ( 1271 ) implies the regularized Faraday equation (fTTl). 
Obviously 

div div(a g) b — b Cg a) = 0. 



Thus equation ( TTTT ) implies that d t div B = and then div B = provided that 
divB|t = o = 0. Notice that due to this formula also the following property holds: 



div (3 = 0. 



(28) 



It remains to derive the regularized energy balance equation f llQh that is the 
most cumbersome point. Following Mill , we write down 



where 



Then 



4(^ + i|B| 2 ) = -( e -g), 

e = div [(E + p+ |B| 2 ) u - (u • B)B] - pu • F. 



(29) 



d t [(E + p+ |B| 2 )u] = [-(e - Q) + d t p + d t B ■ B]u + (E + p + |B| 2 ) d t u. 
By virtue of 01, B and (22) we find 



rd t [(E + p + |B| 2 ) u] = -r {div [(E + p + |B| 2 ) u] — pu • F 
div[(u • B)B] - + • B - Q}u- + |B| 2 ) - T - div(pu)u^) 



= - (E + p+ |B| 2 ) w-r |v 



|B| 



pu — pF ■ u 



div[(u • B)B] -d t p + P-B-Q\u 



Furthermore 
^ E + p+ |B| 2 



= V £ + (Vu)u + Iv (p + ||B| 2 ) + i V (i|B| 2 ) - ^^Vp. 

p p p 



By virtue of ( 1231 ) and ( 1251 ) we also have 

-r4[(u • B)B] = (w ■ B)B + (u • p)B + (u • B)/3. 

Consequently, also subtracting and adding the term r div(B®B) -u, rearranging 
the summands and recalling formula (EH]), we obtain 

r4 [(£■ + p + |B| 2 ) u - (u • B)B] = - (E + p + |B| 2 ) w 



+(w • B)B - rMu - t 



p |^uVe - -^uVp j - Q 



u 



(pw • u)u - r[(u • /3)B + (u • B)/3] + r(-d t p + /3-B)u , (30) 



where 
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M = div(B B) • u - div[(u • B)B] + uV (^IBI 2 ) - S-uVp. 

p 

Differentiating easily leads to formulas 

div(B <g> B) • u = div[(u • B)B] - (BV)u • B, 
uV(i|B| 2 ) = (uV)B-B, 
(3 = (divu)B + (uV)B - (BV)u (31) 

(the last one uses the property divB = 0). Thus 

IBI 2 IBI 2 
M = /3-B- (divu)lBl 2 - l —^-uVp = (3 B - L_ Ldiv(pu). (32) 

P P 

Applying the formula 

(a <S> b)u = (b • u)a, 

from ( 1T51 ) we find that 

Ilu - U NS u = (pw ■ u)u - r[(u • B)/3 + (u • (3)B] 

+r f uVp + pC 2 s divu- — Q + /3-B J u. (33) 
V pe J 

Recalling formula fl23p, we recognize that this term represent one in the curly 
brackets in (1301). 

Therefore inserting formula ( 1321 ) (and recalling (1221)) and ( 1331 ) into (1301). we 
obtain 

rd t [(E + p+ |B| 2 ) u- (u • B)B] = -(E + j9)w - |B| 2 w 
+ (w • B)B + q + xVO - t{(3 ■ B)u - (Ilu - n^u), 

where q is given by formula (TT71). This together with ( 1201 ) straightforwardly 
leads to the regularized energy balance equation (TTOj). □ 

The regularized MGD system of equations (1S|)-(1TB|) generalizes to the mag- 
netic situation the quasi-gas dynamics system in the case of real gas from [7]- 
[10J. It also generalizes the quasi-MGD system from [16] to the case of real gases 
in presence of a body force and a heat source. Notice that our system is written 
in another more standard form that can be essential for further discretization. 
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Proposition 2. For the regularized MGD system of equations, the following 
internal energy balance equation holds 

df(ps) + div[pe(u — w)] + p div(u — w) = div[— q + r(u • B)/3] 
+n : Vu + wVj9 + [- div(B ® B) + V (||B| 2 ) - pF] ■ w 

+r[(uV)B • (3 - (/3V)B -u] + Q. (34) 

Proof. For the regularized MGD system of equations, we subtract the impulse 
balance equation (JSj) multiplied innerly by u and the Faraday equation ( TTTj ) 
multiplied innerly by B from the energy balance equation (HOD. Since 

dt(pu) ■ u + div[p(u — w) <S> u] = dt (|/?|u| 2 ) + div [|/?|u| 2 (u — w)] 

taking into account the mass balance equation (JED, we get 

d t (pe) + div[(pe +p)(u - w) + |B| 2 (u - w) - ((u - w) ■ B)B] 
+ div(B <8) B) • u - uV (p + ||B| 2 ) - div[(u - w) ® B - B <8> (u - w)] • B 
= div [-q + t(/3 ■ B)u + IIu] - div II • u - div[r(u <8> (3 - (3 <8> u)] • B 

-[pw-r div(pu)u] • F + Q. (35) 

The following formulas are valid 

div(v <8> B - B <8> v) • B 
= div [|B| 2 v - (v • B)B] - (vVB) • B + (BV)B • v 
= div [|B| 2 v - (v ■ B)B] - vV (||B| 2 ) + div(B ® B) ■ v 

since div(B ® B) = (BV)B due to divB = and 

div[r(u ® (3 - (3 <8> u)] • B 
= div[r(/3 • B)u - r(u • B)/3] - r[(uV)B • /3 - (/3V)B ■ u]. 

Exploiting them in ( 1351 ) for v = u — w, differentiating in the terms div[j?(u — 
w)] and div (IIu) and applying ( 1221 . we derive 

dt(pe) + div[/?e(u — w)] + pdiv(u — w) — wV]9 
= div[-q + r(u • B)/3] + II : Vu + [- div(B B) + V (||B| 2 )] • w 
+r[(uV)B • (3 - (/3V)B • u] - pw ■ F + Q. 

This equality implies the internal energy balance equation (1541. □ 



S 



Notice that the right-hand side of equation (Elj) does not depend on B for 
t = 0, i.e. without the regularization (and that is the reason why the right-hand 
side of equation ( I24f ) does not depend on B too). 

The crucial point of the physical correctness of a gas dynamics system is an 
adequate entropy balance equation. The entropy s = s(p,e) can be introduced 
by the Gibbs formulas 

p 1 

W Se = e' 

see [IB] • The next proposition generalizes to the magnetic situation the corre- 
sponding results from [71 El ID] • 



S P = s £ = 77> ( 36 ) 



Proposition 3. For the regularized MGD system of equations, the following 
entropy balance equation holds 

d t (ps) + div[ps(u - w)] = div (-|) + - Q S, (37) 
where the entropy production — S is expressed by a formula 

- - , P\~\2 , t Pp\a- ( m2 , rpe e (9p e Q \ 2 
n = n NS ,o + - w H y - div(pu) — divu + uW-- 

+ ,| / 3P + (B.u)^ + 0(l- 3 l|-) (38) 

with -Sjv5,o being the Navier- Stokes entropy production for Q = 0: 

S^ s> o = 2/2ByDy+ ( A "^) ( divu ) 2 + |l V ^! 2 ^° /or n = 1,2, 3. (39) 
Tae entropy production can be also expressed by a formula 

S = + - |w| 2 + ( /9C s 2 divu + uVp - J 

7 7 ^ ^uVs - ^- ) + r\f3\ 2 + (B • u)/3Vr + Q ( 1 - -^-] (40) 



provided that p p > 0. #ere c p and c w = e# are £ae specific heats of gas at 
constant pressure and at constant volume related by a formula [18] 

Cs _ C P 
Pp C-v 
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Under conditions 

tQ 

r = const, — — ^ 1, (41) 

the entropy production is nonnegative: — S ^ for n = 1, 2, 3. 

Proof. Extending the argument of [9] (see also [71 ED]), we introduce the total 
time derivative 

D t ip = d t {pip) + div[p(p(u - w)] = p5 f ^ + p(u - w)Vy>, 

see the mass balance equation (jHj). In a standard manner we write down 



oil 1 
As = 3 A- + t = 3 [pdiv(u - w) + Ae] 

opt) V 



= div 



9 6 



1 



^ + -(B-u)/3 + -E (42) 
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according to the Gibbs formulas ( 15BI ) and the internal energy balance equation 
(1541), where 

S = i [-q + T (B • u)/3] V0 + n : Vu + wVp 

u 

+ [- div(B <8> B) + V (||B| 2 ) - pF] • w + r[(uV)B • /3 - (/3V)B • u] + Q. 
Furthermore, we have 

n : Vu = U NS : Vu + p(uV)u • w + r f u Vp + pC] div u - — Q ) divu 

V ' P £ e J 

+r [-(/3V)u • B - (BV)u • (3 + (/3 ■ B) divu] 

and 

-xW • V0 + U NS : Vu = E NS)0 . 

Next clearly 

p(uV)u • w + wVj9 + [- div(B <g> B) + V (||B| 2 ) - pF] ■ w 

= ^|w| 2 + - div(/?u)uVj9, (43) 
r p 
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see formulas ( |T4"f ) and (J22D- From [91 [TP] the following representation is known: 



1 

T 



p yaVe - uVpj - Q 
+r I uVp + pCg div u — 



r 



uV/9 H — div(/m)uVp 
P 



Pe 
pee 



Q divu + Q 



= — - divfpuj I H — — div u + u\/0 

P \pe 



Q 



+ Q l 



• (44) 



It remains to collect all the other magnetic terms (containing B) with the 
multiplier r. We have 

M = l(B- u)/3V0 - (/3V)u • B - (BV)u • /3 + ((3 ■ B) divu 



+(uV)B • (3 - (/3V)B • u 



1 



= -(B • u)/3V0 + [-(BV)u + (divu)B + (uV)B] • (3 - (/3V)(B • u) 

V 

= i(B-u)/3V^+|/3| 2 -div[(B-u)/3] 



= \(3\ 2 -e&iv 



IB ■ u)/3 

; 



where formula ( 1HTT ) and property ( 1251 ) have been applied. Consequently 



div 



r 



r 



1 



-(B • u)f3 +- Q M = - (r\f3\ 2 + (B • u)/3Vr) 



Formulas 



e K \ o e 

imply the stated ones 



and 



Another representation ( 139| ) for S also follows from [7] , P [10] . 



(45) 



□ 



Notice that formulas ( 1381 ) and ( jlQj ) remain valid in the case r ^ provided 
that one rewrites the term £ Iwl 2 as 



T 



|p(uV)u - div(B ® B) + V (p + ||B| 2 ) - pFp 
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